Tan's distributions and Fermi-Huang pseudopotential in momentum space 
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The long-standing question of finding the momentum representation for the s-wave zero-range in- 
teraction in three spatial dimensions is here solved. This is done by expressing a certain distribution, 
introduced in a formal way by S. Tan [Ann. Phys. 323, 2952 (2008)], explicitly. The resulting form 
of the Fourier-transformed pseudopotential remains very simple. Operator forms for the so-called 
Tan's selectors which, together with Fermi-Huang pseudopotential, largely simplify the derivation of 
Tan's universal relations for the Fermi gas, are here derived and are also very simple. A momentum 
cut-off version of the pseudopotential is also provided, and with this no apparent contradiction with 
the notion of integrals in Tan's methods is left. The equivalence, even at the intermediate step level, 
between the pseudopotential approach and momentum-space renormalization of the bare Dirac delta 
interaction is then shown by using the explicit form of the cut-off pseudopotential. 
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I. INTRODUCTION 

The physics of strongly interacting quantum many- 
body systems has been one of the most active fields of 
research for a number of decades now These sys- 
tems, which are relevant in many areas of physics, such 
as atomic, molecular and optical 0, condensed matter 
3 or nuclear physics f3| , offer challenging theoretical and 
experimental problems that are attracting much of recent 
interest. 

Lately, one of the most studied problems is that of 
a spin-1/2 many-fermion system interacting with short- 
range potentials with a large two-body scattering length, 
i.e. close to unitarity. When this is the case, perturba- 
tion theory is not useful due to the absence of an obvious 
small parameter in the system, and solving the problem 
becomes a major challenge. In this direction, quantum 
Monte Carlo and density-functionals as well as the 
epsilon expansion have been succesfully applied. In 
addition, some of the most striking properties of inter- 
acting Fermi gases at unitarity are those related to its 
universality [7'], due to the absence of a natural length 
scale associated with the interaction. 

Close to but away from unitarity, the so-called Tan's 
relations dictate the behavior of the many-fermion 

system. These important results relate many of the 
many-body properties, such as the adiabatic change of 
the energy when varying the two-body scattering length, 
the asymptotic momentum distribution and the pressure 
to a single quantity called the contact, and were recently 
verified experimentally by Stewart et al. llj. In his pi- 
oneering works, S. Tan introduced, quite formally, the 
regularized s-wave zero-range, or Fermi-Huang [l2l . [l3| 
pseudopotential in momentum representation. He also 
introduced a set of distributions - selectors - in the same 
fashion of his approach to the interaction. These dis- 
tributions proved very useful in the derivation of Tan's 
relations, radically simpler than with conventional ap- 
proaches used later, such as the operator product expan- 



sion [1J|. The tools developed in [8| seem, at first sight, 
to run into contradiction, for a new notion of improper 
integrals or the existence of new types of generalized func- 
tions - distributions - was needed in order to introduce 
them. 

In this Letter, by explicitly constructing Tan's gen- 
eralized functions, we find that they are nothing but 
usual distributions containing Dirac deltas. Moreover, 
related distributions - which become the desired gener- 
alized functions after taking limits - are constructed so 
that integrals over the whole momentum space can be 
defined as limits of cut-off integrals, as in "usual" math- 
ematics. These findings put an end to any kind of con- 
troversy or contradiction regarding Tan's approach [l^ . 
Moreover, they solve the long-standing problem (54 years 
old!) of finding the Fourier-transformed Fermi-Huang 
pseudopotential, which may be useful for few- and many- 
body problems worked out in momentum representation, 
and for variational or perturbative calculations done in 
this way. As an important corollary of these results, it 
is proven that renormalization of the Dirac delta inter- 
action, and Fermi-Huang pseudopotential in momentum 
space are equivalent to one another, even at intermediate 
steps in calculations. 



II. FERMI-HUANG PSEUDOPOTENTIAL 

In order to have a model potential which is simple 
enough, still being able to reproduce the low-energy 
scattering properties in a system, the regularized zero- 
range s-wave (Fermi-Huang) pseudopotential [l^j is use- 
ful. This two-body interaction has the form 



V{r - r') = gS{r - r') 



d 



d\r 



|r-r'|.). 



(1) 



The differential term multiplying the Dirac delta has the 
effect of removing l/|r — r'|-type singularities at the con- 
tact between two particles. Actually, the derivative in 
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(H]) suggests that if it was absorbed into the definition of 
the couphng constant g, then this would depend on the 
relative momentum k (operator [l6l | ) , probably in a very 
special way. This will clearly be one way of interpreting 
the results found in the following section. 



III. PSEUDOPOTENTIAL IN MOMENTUM 
SPACE 

In Tan's pioneering work 8], the formal momentum 
representation of the s-wave Fermi-Huang pseudopoten- 
tial S{r)dr{r-) is obtained by defining the A-distribution 
via 



<5(r)|-(re*"-)EE<5(r)A(k), 



(2) 



where A is defined as A(k) — 1 for k < oo together with 
the following integral relation 



J dkfc-2A(k) = 0. 



(3) 



A naive, rather illuminating interpretation of the above 
equations is that A is equal to 1 everywhere except at 
infinity, where it behaves as a distribution and not as a 
usual function. Moreover, relation (jS]) has a clear physical 
meaning: the two-body scattering length a diverges at 
infinitely large coupling constant g; indeed, as is well 
known, g (x a. 

The choice for A that can be made to satisfy both Eqs. 
@ and ® has the form A(k) = l-g{k)S{T{k)), where G 
and are ordinary real functions of the momentum, and 
have a singularity at k = 0. Introducing this choice into 
(O, we obtain 5(k) = -7^(k) = l/k, where the integral 
in Eq. ([3]) is understood after performing the change of 
variables l/k = t. The explicit form of the A-distribution 
is therefore given by 



A(k) = A{k) = 1 - 



<5(l/fc) 



(4) 



The L-selector is defined as zero for any fc < oo, 
but satisfies 



dk i(k) 



= 1. 



(5) 



Following a similar analysis as for the A-distribution, we 
obtain the explicit form of the L-selector, 



L(k) = 27r 



^6il/k) 



fc2 



(6) 



and all its properties, derived by Tan in [8|, follow easily 
from the above relation. In addition, its position repre- 
sentation £{r) is given by 



^(r)^ 



^L(k)e-r^-4..(ry|. (7) 



The 77-selector, defined as 77(k) = A(k) + L(k)/2/ig, with 
fi the two-body reduced mass, has therefore the following 
form 



r?(k) = 1 



1 



^gk'- 



■(M<?A:-7r2)5(l/fc). 



(8) 



We are now in position of writing down, for instance. 
Tan's energy theorem Q for a (homogeneous) spin- 1/2 
Fermi gas with contact interactions between particles of 
different spin, totally explicitly as 



k,(T 



Sim 



(9) 



with E the energy of the system, the sum in a running 
through the two possible spin components and m — 
2/i the single-particle mass. The Dirac delta in the sum 
has to be regarded, of course, as a limit representation 
before the thermodynamic limit is taken where the sum 
becomes an integral. 



This is the main result of this Letter. In particular, it 
means that the momentum-space representation of the 
two-body pseudopotential has the form V{'k,'k')/g — 
l--(5(l/|k'|)/|k'|; V is obviously non-Hermitian - and so 
it is with A as defined by Tan Q , Eq. ^ - but this repre- 
sents no problem and is analogous to what happens to the 
one-dimensional hard-sphere Bose gas recently discussed 
in [1^. Eq. puts an end to the long-standing quest 
of finding the explicit form of the Fermi-Huang pseu- 
dopotential in momentum space. It is now clear what A 
actually does. First of all, it lets the interaction be the 
bare zero-range two-body potential at finite momenta, 
while the singular, momentum-space delta-like behavior 
removes the ultraviolet divergence of integral ([3]) at in- 
finity. 

Tan also introduced a set of useful distributions Q, 
which he called selectors L{\i) and ?7(k), that fundamen- 
tally simplify the derivation of his celebrated relations 



IV. CUT-OFF MOMENTA AND 
WELL-DEFINED LIMITS OF IMPROPER 
INTEGRALS 

In practice, e.g. in numerical calculations, a momen- 
tum cut-off kc < 00 is often required. However, the pseu- 
dopotential in momentum space needs the evaluation of 
integrals over infinite space, as seen in Eq. ([3])- This fact 
triggered the claim Q that a new notion of integrals was 
needed to define Tan's distributions. This problem is al- 
ready absent if we use Eqs. (0]) and ^ but, in order to 
make finite integrals with the pseudopotential have a sen- 
sible meaning, we need to define a cut-off A-distribution 
in the following way: 



A'''={k) = 9{kc-k) 



Sjl/k-l/k,) 
k 



(10) 
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where 9{k) is the Heaviside step function. Obviously, 

A = hnife^^^oo A*^": . It is easy to see that A'^^ = 1 for 
k < kc and 

[ dkk-^A''-'{k) = [ dkfc-2A'^-(k) = 0, (11) 

which is the cut-off version of (jSj. Similar cut-off ver- 
sions for the L and rj selectors are trivially obtained in 
the same way. From Eq. (fTTj) . we see that the "prob- 
lem" of infinite integrals not being limits of finite inte- 
grals [1] is non-existent if we use A'^'^. Intuitively, this 
cut-off-dependent interaction brings down the infinity to 
the value of the cut-off chosen for a particular compu- 
tation. It may also be replaced by a convenient limit 
representation with usual functions in numerical calcula- 
tions. 



EQUIVALENCE WITH MOMENTUM-SPACE 
RENORMALIZATION 



The first of the integrals in the second line of Eq. p4)) 
is the one appearing in the ill-defined problem of find- 
ing the bound state for the (irregular) delta interaction 
.9s<5(r), 



= lim 



dk 



k^/2^i+\E\ 



(16) 



lim 



kc — y^2fi\E\ arctan 



fee 



Here, the interaction strength plays the role of the so- 
called bare coupling constant. The rcnormalized coupling 
constant gn is then defined as 



1 



1 



(17) 



and the limit fee — > oo is taken while keeping gn constant. 
We obtain a cancelation of the linearly divergent term kc, 
and finally 



We proceed now to see how the above procedure is 
totally equivalent to the usual renormalization scheme 
for the delta interaction. Rather complete treatments of 
renormalization with effective interactions are given in 
(l"7l - [l9| . Consider the stationary Schrodinger equation 
H-tjj — El}) in the relative coordinate, with 



p2 d 

H='^+g5{r)-{v). 
zji or 



(12) 



We consider bound-state {E < 0) solutions in the mo- 
mentum representation, and we set h = 1 throughout. 
The integral equation for the energy E — —\E\ has the 
form 



1 
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dfc- 



fc2A(fc) 
'k^/2fj,+ \E\ 



(13) 



In order to evaluate the above integral without making 
use of the formal properties of A Q, we note that, using 
Eq. (UnD, 



dfc 



PA{k) 
k^/2^i + \E\ 



(14) 



lim 



dfc 



fc2 



ky2^i+\E\ 



S{t - \/kc)T- 



\E\t^ 



where we have changed variables as fc = l/r in the second 
integral. From Eqs. (|T3l) and ([M]) we obtain 



lim 

- ^ 



^2ijl\E\ arctan 



- fcn 1 - 



1 



k? 



(15) 



which gives, after using g = 27ra//x, the bound state en- 
ergy E = —l/2fia^, as should be. 



E : 



2n^ 



(18) 



It is now obvious that if the rcnormalized coupling con- 
stant gn coincides with the original Fermi-Huang cou- 
pling constant g, then Eq. ([18]) becomes E = —1/2/10^, 
and coincides with the energy calculated with the A- 
distribution. Actually, from Eqs. (|15m7|) . we observe 
that the use of Fermi-Huang pseudopotential in momen- 
tum space and the renormalization procedure defined in 
([T7| are equivalent - except for a negligible factor of 
(1 -|- 2fi\E\/kc)~^ which we could arbitrarily include in 
the renormalization - provided that gu, = 9- This last 
fact is of course natural, since the bound state energy 
E — —\/2^a? is the fit parameter used in renormal- 
ization approaches; actually, both renormalization and 
Fermi-Huang pseudopotentials aim at fitting the bound- 
state energy (or equivalently the scattering length). 



VI. PERTURBATION THEORY 

We show now how perturbation theory can be imple- 
mented using the A-distribution without the need of per- 
turbative renormalization, and that such two approaches 
are again equivalent to one another. 

Let us consider two-body scattering from the Fermi- 
Huang pseudopotential, described in position represen- 
tation by Hamiltonian (1121) in the relative coordinate. 
The system's t-matrix T{z), z e C, is the solution to the 
Lippmann-Schwinger equation T[z) — V + VG''^^T{z), 
withG(°) = G'(°)(z)the non-interacting Green's function. 
In momentum space, the matrix elements (p'| T{z) |p) = 
T(z;p',p) we calculate are at finite momenta p' and p. 
Therefore the only elements that require the regulariz- 
ing term S{l/k)/k are those which involve fc oo, that 
is, only (p'|y|k) = V(p','k) when it is in an integral 
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over dk. Taking this fact into account, we can calculate 
the perturbative expansion (Born series) for the t-matrix, 
with all its terms being finite. Note that due to the spe- 
cial form of the contact interaction, T(z; p', p) = T{z) is 
independent of p', p. 

On-shell (z = 9^/2/^ + it]), the first order term t*-^-* 
of the Born series is given by t^^^ = ^(p',p) = 9- It 
is identical to the first Born obtained for a bare delta 
interaction after first-order renormalization gs ^ gR — g 
(which involves no infinities). The second order term t'^' 
is easily calculated as 



dkk^ 



1 - k-^S{l/k) 
{q^ - k^)/2fi + ir]' 



(19) 

The above quantity is found to be finite, and is given 
by ifj,g^q/2TT^ . The reader can easily check that second- 
order perturbative renormalization for a bare Dirac delta 
interaction yields exactly the same result, provided that 
the renormalized coupling constant gji = g. We may 
continue the perturbative expansion to arbitrary order 
to find that all terms t("\q^ /2n + irj) = {i^igq/2TT)'^-^g 
coincide with those obtained through perturbative renor- 
malization, and resum the resulting series to obtain the 
exact t-matrix T{q^ /2ii + irj)/ g = (1 — i^gq/2TT'^)~^ . 

There is yet another way of calculating the t-matrix 
perturbatively. First note, from its exact form, that 
it is holomorphic in the inverse coupling constant A = 
{2iT)^/g. We can therefore expand it in powers of A. We 
may also obtain the series by calculating exactly only 
the t-matrix at infinite coupling constant T^{z). This 
may be as well considered as a different renormalization 
prescription, where the on-shell Too is measured in an ex- 
periment, although regularization- renormalization of in- 
finities is here not required. Since A = V/g^ we have 
AG"roo(z) = -1, which is solved on-sheh by T^{q^ /2ijl+ 
irj) = 2n^i/fiq; formally, T^{z) = -{AG^°^)-\ We im- 
mediately obtain T — Too — ATooT, which is easily solved 
perturbatively by T = Too X]n(~-^^oo)"- After resum- 
mation the t-matrix is given by T = Too/ (1 + ATqo), the 
exact result also obtained above. 



VII. CONCLUDING REMARKS 

We have obtained the explicit momentum-space rep- 
resentation of the zero-range s-wave pseudopotential and 
Tan's selectors. Their forms are very simple and involve 
usual Dirac delta distributions only. This implies that 
the generalized functions defined by Tan in [8| do not 
constitute a novel kind of distributions, therefore clearly 
justifying the methods he used to derive his celebrated 
relations for the Fermi gas. We have also con- 

structed a cut-off version of the pseudopotential, which 
may be useful for finite-size computations in momentum 
space. The cut-off pseudopotential is crucial to make 
improper integrals involving Tan's distributions become 
limits of integrals over a finite interval, further justify- 
ing Tan's methods. We have then shown that the use of 
this potential or a renormalization procedure are totally 
equivalent, both approaches aiming at removing ultravi- 
olet divergences, and we have exemplified this with the 
two-body bound state problem and perturbation theory. 

Our results constitute a promising path towards ob- 
taining the momentum representation of higher partial 
wave pseudopotentials [l^, which can be of great inter- 
est for the many-body problem in momentum space with 
such interactions; this could simplify many calculations, 
as in the presently studied case of s-wave interaction. An 
important implication of the results of this Letter is that 
the derivation of many celebrated results in the literature 
may become more accessible for a non-specialized audi- 
ence having only a working knowledge of or a preference 
for explicit operator methods. 

It would be very interesting to investigate the explicit 
form of three-body contact interactions for bosons in 
momentum space leading to the equivalents of different 
three-body renormalization schemes [2l|, . 
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